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ABSTRACT 

Baryonic acoustic oscillations imprinted in the galaxy power spectrum provide a promising tool for 
probing the cosmological distance scale and dark energy. We present results from a suite of cosmological 
N-body simulations aimed at investigating possible systematic errors in the recovery of cosmological 
distances. We show the robustness of baryonic peaks against nonlinearity, redshift distortions, and mild 
biases within the linear and quasilinear region at various redshifts. While mildly biased tracers follow 
the matter power spectrum well, redshift distortions do partially obscure baryonic features in redshift 
space compared to real space. We calculate the statistical constraints on cosmological distortions from 
N-body results and compare these to the analytic results from a Fisher matrix formalism. We conclude 
that the angular diameter distance will be constrained as well as our previous Fisher matrix calculations 
while the Hubble parameter will be less constrained because of nonlinear redshift distortions. 



Subject headings: cosmological parameters ~ 
distance scale — methods: 



- large-scale structure of universe 
N-body simulations 



cosmology: theory 
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1. INTRODUCTION 

Baryons create a distinct oscillatory signature in the 
power spectrum of the large-scale structure of the uni- 
verse (Peebles & Yu 1970; Bond & Efstathiou 1984; Holtz- 
man 1989; Hu & Sugiyama 1996; Eisenstein & Hu 1998a). 
These baryonic acoustic oscillations have been seen in the 
anisotropics of cosmic microwave background (Miller et 
al. 1999; de Bernardis et al. 2000; Hanany et al. 2000; 
Halverson et al. 2002; Lee et al. 2001; Netterfield et al. 
2002; Benoit et al. 2003; Bennett et al. 2003; Pearson ct al. 
2003) and recently in large galaxy redshift surveys (Eisen- 
stein et al. 2005; Cole et al. 2005). Because the density 
of baryons is less than that of cold dark matter, the oscil- 
lations in the matter power spectrum are weaker in am- 
plitude than those in the cosmic microwave background 
(hereafter CMB). However, whether in the CMB or in late- 
time structure, the oscillations define a constant comoving 
length scale in linear perturbation theory. 

In our previous paper (Seo & Eisenstein 2003, hereater 
SE03), we demonstrated that large galaxy redshift surveys 
can constrain the Hubble parameter and angular diame- 
ter distance to a precision of a few percent using the im- 
printed baryonic acoustic oscillations as a standard ruler. 
The physical scale of the oscillations can be determined 
from matter density and baryon density of the universe, 
which in turn are deduced from the shape and relative 
amplitude of the baryonic peaks in CMB anisotropy data 
(Eisenstein et al. 1998; Eisenstein 2003). One can then 
compare the physical scale with the observed length scales 
of oscillations in transverse and line-of-sight directions in 
galaxy redshift surveys to yield the angular diameter dis- 
tance and Hubble parameter at the given redshift. This 
in turn measures the evolution of dark energy as well as 
the spatial curvature and Vim- See Blake & Glazebrook 
(2003), Linder (2003), Hu & Haiman (2003), Amendola et 
al. (2005), Cooray (2004), Dolney et al. (2004), and Mat- 
subara (2004) for similar studies. 

Because of its weak amplitude modulations in matter 



distribution, the baryonic features are susceptible to the 
erasure by the nonlinear coupling of Fourier modes (Jain 
& Bertschinger 1994; Meiksin, White, & Peacock 1999; 
Meiksin & White 1999; Scoccimarro et al. 1999), such as 
can result from nonlinear growth, nonlinear bias (Kaiser 
1987; Coles 1993; Scherrer & Weinberg 1998; Dekel & La- 
hav 1999; Coles et al. 1999; Seljak 2000; White 2001) or 
nonlinear redshift distortions (Hamilton 1998; Hatton & 
Cole 1998; Seljak 2001; White 2001; Scoccimarro 2004). 
Nonlinear mode-coupling also produces additional small- 
scale power, thereby altering the shape of the power spec- 
trum and potentially disguising the locations of baryonic 
features. It is therefore crucial to examine whether the os- 
cillatory features are robust against various nonlinearities 
and whether the remaining ripples can effectively distin- 
guish a dilation in distance from other non-cosmological 
effects in the clustering of galaxies. 

In this paper, we present an N-body study of the effect 
of nonlinear growth, nonlinear redshift distortions, and 
halo bias on the detectability of acoustic oscillations. Pre- 
viously, Meiksin, White, & Peacock (1999) used N-body 
simulations to show the effects of nonlinearity on baryonic 
signatures in the present-day large-scale structure. They 
studied the effect of bias and redshift distortions for var- 
ious cosmological models. Our study extends their work 
to higher redshifts. We adopt a ACDM model consistent 
with Wilkinson Microwave Anisotropy Probe (WMAP) 
data (Spergel et al. 2003) and generate density fields at 
redshifts of 3, 1 and 0.3. We then investigate and quantify 
the erasure of baryonic features in the matter power spec- 
trum and biased galaxy power spectrum in real space and 
redshift space at those redshifts. We also attempt to re- 
move the nonlinear alteration from the shape of the power 
spectrum and thereby recover the underlying contrast of 
baryonic features. As this work was being completed, re- 
lated studies by Springel et al. (2005) and Angulo et al. 
(2005) appeared on this topic. 

In SE03, we used the Fisher information matrix to cal- 
culate predictions for the statistical constraints on dark 
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energy. We used a conservative choice of nonlinear scale 
^max {=tt/2R) by requiring an ~ 0.5 and excluded power 
in smaller scales from our analysis. For larger scales, we 
adopted a linear growth function, linear redshift distor- 
tions, and a linear bias model with an additive offset. In 
reality, as the transition from linear to nonlinear scales is 
not discrete, the effects of nonlinearity, i.e., mode coupling 
and more complicated scale-dependence, may mildly con- 
taminate the power spectrum even on large scales. Here, 
wc use our N-body simulations to assess the impacts of 
these nonlinear effects on the baryonic features on large 
scales relative to the statistical constraints we calculated 
in SE03. We compare our N-body results with the choices 
of nonlinear scale in SE03. Our results will provide further 
guidance for linear approximations in various studies with 
baryonic physics. 

Finally we want to investigate the distance constraints 
available from galaxy surveys taking account of the full N- 
body effect. The nonlinear effects not only inhibit us from 
detecting the weak baryonic signatures on small scales but 
also increase the statistical variance of power spectrum 
above the Gaussian estimates. We perform analysis on 
our simulated power spectra to fit to the cosmological dis- 
tances and compare the constraints with those in SE03. 
The result will show whether we can deduce the informa- 
tion on cosmological distances from the power spectrum 
altered by the nonlinear effect. 

In § 2 wc describe the parameters of our cosmological 
N-body simulations. In § 3 we present the effect of nonlin- 
ear growth and redshift distortions on baryonic features in 
the matter power spectrum. In § 3.3 we remove the effect 
of nonlinearity on the broadband shape and study the re- 
sulting features. In § 4 we study the baryonic signatures 
in biased power spectra in real space and redshift space. 
In § 5 we present the errors on cosmological distances re- 
sulting from analysis. 

2. COSMOLOGICAL N-BODY SIMULATIONS 

We run a series of cosmological N-body simulations us- 
ing the Hydra code (Couchman, Thomas, & Pearce 1995) 
in colhsionless P^M mode. We use the CMBfast (Seljak & 
Zaldarriaga 1996; Zaldarriaga et al. 1998; Zaldarriaga & 
Seljak 2000) linear power spectrum to generate many ini- 
tial Gaussian random density fields at redshift of 49 and 
evolve them to redshifts of 3, 1, and 0.3. The cosmological 
parameters wc use to generate the initial fields are flm = 
0.27, nx = 0.73, Oh = 0.046, h = 0.72, and n = 0.99. The 
initial fields are normalized by requiring cts = 0.9 at z = 
and assuming a linear growth function. Each simulation 
box represents T^ox = 512^h~^ Mpc^ and contains 256^ 
dark matter particles 8.28 x lO^^Afsun/particlc). We 
compute gravity using 256"^ force grids with a Plummer 
softening length of 0.2 Mpc. We use 51 simulations 
at .2 = 1 and 0.3, and 30 simulations at z = 3. Note that 
the total volume of the simulations is much larger than 
the survey volume parameters listed in SE03 so that we 
can study the effect of nonlinearity with little interference 
from statistical variance. 

The resulting density field of each simulation box is 
Fourier transformed, and the squared complex norms of 
Fourier coefficients are spherically averaged over all simu- 
lations to give the matter power spectra in the wavenumber- 
shells of widths Ak = 0.005/i Mpc~^. A mode is in- 



cluded in a shell if its discrete wavenumber falls in the 
shell. Because each mode in the discrete transform is in- 
cluded in one and only one wavenumber-shell, the shells 
are not correlated for a Gaussian field even though the Ak 
is smaller than the size of an independent cell in Fourier 
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space, ^tt/V^^^^ . Using small AA; ensures that the contrast 
of narrow features are not artificially reduced. However, 

because a thin shell contains fewer modes, wc do want to 
apply some smoothing. We use Savitzky-Golay filtering 
(Press et al. 1992), as this can preserve peak heights bet- 
ter than boxcar smoothing. The Savitzky-Golay method 
also gives estimates of the derivatives"*^. Prom compar- 
isons between power spectra before and after smoothing, 
we believe that this procedure neither introduce mislead- 
ing baryonic features nor erase meaningful features. For 
power spectra in redshift space, we spherically average the 
Fourier transform of the density fields after displacing par- 
ticles according to their peculiar velocities assuming a dis- 
tant observer. 

3. THE NONLINEAR MATTER POWER SPECTRUM 

3.1. Nonlinear effects in the matter power spectrum 

As mass perturbations on a given scale approach or- 
der unity in amplitude, linear perturbation theory breaks 
down and the gravitational growth of perturbations in one 
mode is increasingly coupled with perturbations in other 
modes. The higher-order contribution resulting from this 
mode coupling hinders the detection of features in the ini- 
tial power spectrum, including the baryonic acoustic os- 
cillations. That is, the additional power contributed from 
other modes blurs the initial features at a given Fourier 
mode as they mix with a convolution of other modes (Jain 
& Bertschinger 1994, and references therein). Nonlinear 
growth from mode coupling increases power above the lin- 
ear growth rate on small scales, resulting in a bigger sta- 
tistical variance for any underlying initial features. As the 
amplitude of density perturbations grows with time, these 
nonlinear effects become stronger and proceed to larger 
scales. 

As a basic model, one might distinguish the linear 
regime from nonlinear regime with a scale R (= 7r/2A;niax) 
for an appropriate rms overdensity fiuctuation, an, at a 
given epoch, and assume scales with a smaller rms over- 
density as linear. In studies of the statistical expectations 
from large redshift surveys, imposing more conservative 
values of R improves the linear approximation, but this is 
at the expense of more of the remnant linear information 
in the quasilinear regime beyond R. On the other hand, 
more forgiving criteria for R will result in an overestima- 
tion of the performance: the precision of the acoustic scale 
measurement improves as fcmax; increases, finally saturat- 
ing beyond about 0.25/i Mpc~^ because of Silk damping 
(see Figure 4 of SE03). Only an N-body study can say 
whether this nonlinear scale accurately accounts for the 
erasure of features in the initial power spectrum such as 
baryonic oscillations. This will give an additional handle 
in the transition from linear to nonlinear regime beyond 
that traced by the increased amplitude. We will use the 
notation A;niax in this paper to characterize the nonlinear 

^We use 4*'' order polynomial with a filter width of Ak = 
0.04h Mpc~^ to smooth power spectra and Afe = 0.05/i Mpc~^ to 
derive the derivatives of power spectra. 
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Fig. 1. — The matter power spectra at various redshifts. Left: each power spectrum is divided by a zero-baryon power spectrum at the given 
redshift. SoUd Hues are for the real-space clustering, and dashed lines with the same color are for the corresponding redshift-space clustering. 
Right: dlnP/dlnk from the matter power spectra at various redshifts in real space (upper panel) and in redshift space (lower panel). Green 
dashed line: the input power spectrum from the CMBfast, black line: the linear matter power spectrum at ^ = 49 generated from the input 
power spectrum, red: the nonlinear matter power spectra at 2 = 3, blue: 2 = 1 and violet: z = 0.3. The first troughs in dlnP/dlnk from the 
N-body results are lower than that of the input power spectrum due to the interaction of the Savitzky-Golay smoothing with the boundary 
at fc ~ 0. 



scale and to assess the linear information from the bary- 
onic features surviving nonlinearity. 

Figure 1 shows the effect of nonlinear gravitational 
growth on baryonic acoustic oscillations in matter power 
spectra from N-body simulations. The solid lines in the 
left panel of the figure show the spherically averaged real- 
space power spectra divided by a zero-baryon power spec- 
trum^ (Eisenstein & Hu 1998a) at various redshifts. The 
growth rate calculated from rms overdensity fluctuations 
at 16h~^ Mpc, CiQfi-i Mpci is consistent with the linear 
growth rate, while the values of crg^-i Mpc show nonlinear 
effects at lower redshift. The vertical lines denote the non- 
hnear scales, fcmax, adopted in SE03 to satisfy an ^ 0.5 
(fcmax = 0.11/iMpc-i at z = 0.3, 0.19ft. Mpc-i at z = 1 
and 0.53/1 Mpc^^ at z = 3). 

As expected, nonlinear structure formation increases 
small-scale power and obscures small-scale baryonic fea- 
tures. The effect proceeds to larger scales with time. The 
difference in slope between the linear power spectrum at 
z = 49 and the power spectra at lower redshift shows that 
remnant nonlinearity exists even for k < k^a,x- At lower 
redshifts, the contrast of baryonic peaks on large scales is 
decreased because of the nonlinear mode-coupling. 

The right panel of Figure 1 shows logarithmic derivatives 
of the matter power spectrum with respect to wavenum- 
bers, generated from Savitzky-Golay filtering. This plot 
of derivatives is useful not only because it effectively 
manifests fine details of power spectra but also because 
d In P/d In k is what enters into the Fisher information ma- 
trix and creates the standard ruler test with baryonic os- 
cillations. Of course, the derivative din P/ dink increases 
the noise in the power spectra despite our use of smooth- 

power spectrum generated from the fitting formula in Eisen- 
stein & Hu (1998a) for our fiducial Qm and h but fi;, = 



ing. 

The real-space d In P/d In k in Figure 1 (upper panel) 
demonstrates that the oscillatory features are well distin- 
guished for k < fcmax at z = 0.3 and z = 1 despite the 
slight nonlinearity. At z = 3, it becomes difficult to distin- 
guish baryonic features beyond k 0.4/iMpc~^ because 
the decreased contrasts of small-scale peaks caused by Silk 
damping (Silk 1968) make them susceptible to even small 
amount of nonlinearity or noise. Fortunately, baryonic 
features in A: > 0.3ft. Mpc~^ have minor contributions to 
cosmological information (SE03). 

3.2. The effect of redshift distortions 

Since we measure redshifts of galaxies rather than their 
physical distances, the three-dimensional galaxy power 
spectrum is subject to redshift distortions, the angle- 
dependent distortion in power spectra caused by the pe- 
culiar velocity of galaxies (Kaiser 1987; Hamilton 1998; 
Scoccimarro 2004). On large scales, the bulk motions 
of large-scale structure toward overdense regions enhance 
power, and on small scales, the virial motions within and 
among halos create an apparent extension along the line of 
sight, known as the finger-of-God effect (de Lapparent et 
al. 1986). This suppresses power on small scales. In linear 
theory, the large-scale power enhancement by the redshift 
distortions follows a simple form (Kaiser 1987). This is 
true only for the asymptotic limit of large scale, and in 
general, the nonlinear effect in velocity fields deviates the 
redshift-space power from Kaiser formula even on fairly 
large scales (Scoccimarro 2004). We seek to estimate how 
much the nonlinear effect of large-scale redshift distortions 
affects the baryonic features. 

The dashed lines in the left panel of Figure 1 show 
the redshift-space power spectra divided by a zero-baryon 
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Fig. 2. — Left; dhi P/dliik from Pnonlinear — fh{co,k,k^) at various redshifts in real space. Gray solid line: the linear matter power 
spectrum at 2 = 49, black solid lines: the nonlinear matter power spectra at various redshifts. Right: the redshift-space dhi P/dhik from 
-Pnonlinear ~ /L(c(),fc,fe^) after Corrected for the finger-of-God suppression. Note that the gray lines in this panel denote the linear matter 
power spectrum in real space at z = 49 while black lines are for the nonlinear matter power spectra in redshift space. The vertical dashed lines 
denote the value of fcmax we assumed in SE03. We conclude that the baryonic features survive on large scales despite the nonlineaj: growth 
with the redshift distortions imposing an additional but mild degradation. 



power spectrum. The figure depicts tlie progress of redshift 
distortions with nonhncarity and their effect on the bary- 
onic features in the matter power spectrum. The redshift- 
space power spectra on large scales have a higher ampli- 
tude than that of real space by the amount predicted by 
linear theory (Kaiser 1987) for the asymptotic limit. Even 
on large scales, we observe a slight suppression in redshift- 
space power with respect to linear theory, in agreement 
with Scoccimarro (2004). On small scales, the finger- 
of-God effect not only suppresses the overall power but 
also decreases the contrast in baryonic features. As ex- 
pected, the finger-of-God effect increases with time, and 
the resulting suppression makes the position of the second 
peak (which is beyond the nonlinear scale for z = 0.3) 
appear slightly shifted. For our cases of study, the finger- 
of-God effect generally produces 10 — 20% of suppression 
in redshift-space power at fcmax in all redshifts when com- 
pared to the prediction by linear theory. 

The lower right panel of Figure 1 gives a more clear view 
of the redshift distortions decreasing the contrasts and in- 
troducing noise at A; < fcmax- Nevertheless, the baryonic 
features are still preserved up to fc ~ fcmax for z = 0.3 
and z = 1, and fc < 0.3ft- Mpc"-'^ for z = 3. The result at 
z = 0.3 seems especially encouraging in that the features 
are preserved even beyond fcmax- 

It is important to note that these curves are spherically 
averaged power spectra in redshift space while we aim to 
use anisotropic information of power spectra in real ob- 
servations. In the three-dimensional power spectra in red- 
shift space, wavcvcctors nearly perpendicular to the line- 
of-sight direction will preserve baryonic features as well as 
the real-space power in Figure 1, and wavevectors nearly 
along the line-of-sight direction will appear more smeared 
than the redshift-space power in the figure. 



3.3. Acoustic features after restoration of the broadband 

shape 

In this section, we consider the effect of nonlinear mod- 
ification of the slope of the power spectrum on the bary- 
onic features. While the nonlinear mode-coupling effect 
directly erases the baryonic features, the accompanied in- 
crease in small-scale power due to nonlinear growth will 
modify the slope of the power spectrum as a function of 
wavenumber. The resulting change in slope will shift the 
apparent locations of the baryonic peaks. A careful look at 
Figure 1 reveals that the higher harmonics of the baryonic 
peaks appear at slightly larger wavenumbers compared to 
the initial power spectrum. Furthermore this addition of 
broadband power can misleadingly decrease the contrast of 
the peaks and thereby overestimate the loss in information 
by the mode coupling effect on the baryonic features. 

We attempt to correct for the nonlinear growth effect 
on the broadband shape of the power spectrum by mod- 
eling the gradual modification to the broadband shape 
due to nonlinearity as a smooth function of wavenum- 
ber and subtracting off this smooth function from the 
nonlinear power spectrum to restore the oscillatory por- 
tion to its original slope. We fit the nonlinear power 
spectra in real space to a multiple of the linear power 
spectrum, g^Punear, and a 2"'^-order polynomial function 
/l = C() + Ci k + C2 (with and the Ci as constants) 
and then subtract the smooth function /l from the non- 
linear power spectra before calculating dlnP/dlnfc. This 
process increases the fractional variation in P by decreas- 
ing the overall amplitude. In other words, noting that 
dlnP/dlnfc = (1/P) x (rfP/dlnfc), we have decreased P 
but not substantially changed its derivative. Increasing 
the degree of the function /l up to 3''''-order does not 
have a sizable effect on rflnP/rflnfc on large scales. As 
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required, the resulting function /l is nearly constant on 
large scales in all cases. The returned values of are not 
consistent with linear growth rate at lower redshift because 
the fitting process tends to decrease to account for the 
erasure of baryonic features. Since we do not want to over- 
amplify the baryonic features by oversubtracting beyond a 
true nonlinear power, as a sanity check, we calculate and 
subtract /l by fixing to the analytic linear growth rate. 
This reduces the contrast of resulting baryonic features in 
rflnP/dlnfc at z = 0.3 but only by a small amount. 

The left panel of Figure 2 shows the resulting nonlinear 
In P/rf In fc in real space (black lines) in comparison to the 
linear power spectrum at z = 49 (gray lines). The agree- 
ment between the power at z = 3 and the linear power 
spectrum is excellent. We also recover larger contrasts 
in ^; = 0.3 and z = \ cases after the broadband shape 
is restored. Our result shows that it is possible to trace 
baryonic features up to k ^ 0.3h Mpc"'^ at z = 1 and 
k ~ 0.2h Mpc^^ at z = 0.3. The results suggest that the 
performance of baryonic features as a standard ruler will 
be diminished due to the decreased contrast at fc < fcmax, 
but the features preserved beyond fcmax will tend to com- 
pensate for the reduction. 

In short, we find that correcting for the nonlinear effects 
on the broadband shape with a smooth function helps us 
rescue baryonic features by some degree, implying that 
the nonlinear effects are relatively smooth in power. We 
interpret the aniomit we cannot recover as having been 
lost to mode coupling effects. 

For the redshift-space power spectrum, we first correct 
for the finger-of-God suppression. We fit the nonlinear 
power spectrum in redshift space to the counterpart in 
real-space by a multiplicative smooth function in the form 
of Ffog = l/(fc'"cr'" + 1)1/™ where a and m are fitting 
parameters. We find m ~ 2.3 — 2.95. After dividing by 
Pfog we calculate and subtract /l to match the restored 
redshift-space power spectrum to the linear power spec- 
trum. Although Ffog may not necessarily match the con- 
ventional form of an exponential finger-of-God suppres- 
sion, the function behaves properly at the limits of small 
and large wavenumbcr and is capable of characterizing the 
difference between the real-space and redshift-space power 
from our N-body results. The right panel of Figure 2 shows 
the resulting din P/ dink in redshift space (black lines) 
in comparison to the linear power spectrum at z=49 in 
real space (gray lines). The restored rcxlshift-space power 
spectra show slightly larger contrasts in baryonic features 
at lower redshift relative to the uncorrected ones (lower 
right panel of Figure 1). From the comparison between 
the real-space and redshift-space din P/ dink (Figure 2), 
the redshift-space power reasonably traces the details in 
the real-space power on linear and quasilinear scales al- 
though the contrast of the baryonic features degrades in 
redshift space. 

In summary, the real-space (redshift-space) matter power 
spectrum traces baryonic features up to fc ~ 0.4 (0.3)/i Mpc" 
at 2: = 3. At z = 1, /Cmax of 0.19/1 Mpc~i seems a rea- 
sonable choice for the linear approximation in real space 
considering the features preserved on even smaller scales, 
but a lenient standard in redshift space. At z = 0.3, 
kmax ~ 0.11ft. Mpc~i is conservative both in real space 
and redshift space. Adopting a slightly larger fcmax ap- 



pears justified at this redshift. 

4. THE EFFECT OF BIAS 

4.1. Anomalous power 

In galaxy redshift surveys, we do not directly observe the 
real-space matter power spectrum but instead observe the 
distribution of biased tracers of matter in redshift space. 
The assumptions of local bias and Gaussian statistics for 
the density fields lead to a scale-independent bias on large 
scales for the correlation function (Coles 1993; Scherrer & 
Weinberg 1998; Meiksin, White, & Peacock 1999; Coles ct 
al. 1999). Any excess small-scale correlation from biasing 
will appear as an additional constant term in the biased 
power spectrum on large scales. This holds even when 
the matter density field is nonlinear (Scherrer & Weinberg 
1998; Coles et al. 1999; Seljak 2000). The bias on large 
scales is thus scale-independent up to this constant term. 
On smaller scales, the bias will generically deviate from the 
simple approximation of scale independence. We use the 
term 'nonlinear bias effect' in this paper to designate any 
deviation from a simple multiplicative bias in the power 
spectrum. 

As the biased tracers such as galaxies are rare objects 
compared to the underlying matter, the biased power spec- 
trum is subject to a shot noise. Conventionally, one writes 
the shot noise as a white noise {P ~ constant) equal to 
the inverse of the number density of particles. However, 
the effects of the limited number of particles can be more 
complicated than a simple Poisson noise, such as in bi- 
asing schemes where halos are not equally weighted or, 
more enigmatically, the discreteness effects in the process 
of finding halos. For example, the fact that a halo finder 
cannot identify two halos arbitrarily close together means 
that the shot noise will not be white. Therefore it is not 
straightforward to identify the shot noise with the con- 
ventional white noise based on an inverse number density. 
Rather than singling out the shot noise term, we hereafter 
group this term with the nonlinear bias term without dis- 
tinguishing one from the other. Both terms are easily ap- 
proximated as an additive constant on large scales. Here- 
after, 'anomalous power' is used to refer to a combination 
of these two terms. We find that the fractional level of the 
anomalous power compared to the linearly biased power is 
important to track. We hereafter will describe this frac- 
tional level as 'larger' or 'smaller' anomalous power. 

The anomalous power will contribute additional power 
above the linearly biased power spectrum, therefore in- 
creasing the statistical variance of the underlying features. 
The nonlinear bias effect in the anomalous power may 
also induce a mixture of information from different Fourier 
modes and so erase features. 

4.2. Bias schemes 

In generating biased tracers of the matter, we do not 
attempt to reproduce realistic galaxy populations but in- 
stead to refiect an interesting range of models of galaxy 
populations using simple deterministic halo-based biasing 
schemes. The current halo occupation distribution (HOD) 
models suggest a single galaxy per halo for low-mass halos 
above a mass threshold and an additional power-law mean 
occupation for more massive halos while the details vary 
in different studies of galaxy populations (Berlind et al. 
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Fig. 3. — Biased power spectra at 2 = 1 divided by a zero-baryon power spectrum. The left panel is for the number-weighted cases (NUM), 
and the right panel is for the mass-weighted cases (MASS). The label 'm' denotes the minimum group multiplicity for selecting halos. For 
the NUM cases, b ^ 1.7, 2, and 2.5, and for the MASS cases, b ^ 2.4, 2.7, and 3.1 as m increases. The solid lines are for the real-space power, 
and the dashed lines are for the redshift-space power. 



2003; Kravtsov et al. 2004; Zehavi et al. 2004). We reduce 
the complexity by decomposing the HOD models into the 
two extreme bounding cases: one in which halos above a 
mass threshold host a single galaxy and one in which ha- 
los above a mass threshold have mean occupation as linear 
to the halo mass (i.e., the power-law index of unity). To 
create different amplitudes of bias, we apply various min- 
imum group multiplicities (i.e., halo mass-thresholds) for 
both cases. Various superpositions of these trial cases then 
can comprise more complex HOD models. Note that the 
resulting populations themselves already represent super- 
positions of different-mass halos. Therefore our biasing 
schemes will be sufficient for examining the robustness of 
the baryonic features in various galaxy populations be- 
cause, if baryonic oscillations are found to survive in both 
of the extreme cases, it seems unlike that a mixture would 
fair worse. This will be further justified once our results 
can show that different bias models indeed extrapolate to 
linear biasing on large scales. 

We use the friends-of-friends method (Davis et al. 1985) 
and identify halos by adopting a linking length of 0.6 
Mpc and minimum group multiplicities of 4, 10, or 
30 particles. We assign zero galaxy density for the re- 
gions below this threshold. For the regions identified as 
halos with a given threshold, we assign galaxies with two 
schemes: 

1. Number Weighted (NUM), where we assign one 
galaxy per halo, using the center-of-mass position and ve- 
locity. In this case, the number of galaxies does not follow 
the mass of halos, and the information on the virialized 
motions within each halo (the finger-of-God effect) is lost. 

2 Mass Weighted (MASS), where we retain the velocity 
and density structure of the halo by assigning one galaxy 
per simulation particle. This way, the number of galaxies 
follows the mass of the halo, and the finger-of-God effect 
is preserved. 



The minimum group multiplicity, to, of 4 and 10 are 
small compared to 20 — 30 particles usually desired for 
robust halo identification (Somerville et al. 2000; Jenk- 
ins et al. 2001). However, the low multiplicity halos are 
still tracing overdense regions, albeit in a more stochas- 
tic manner. This randomness would not be favorable for 
recovering acoustic oscillations, implying that our results 
are conservative. 

At lower redshifts, both methods are able to provide us 
a sufficient range of bias values of our interest. However, 
at z = 3, the mass resolution of our simulation is too low 
to find halos with masses low enough to yield mild values 
of biases (i.e. 6 ~ 3). For this reason, we use an addi- 
tional bias scheme at this redshift, where the density of the 
tracers approximately follows the matter density squared 
(hereafter, Rho^). In this biasing scheme, each particle is 
weighted by the average density with a 2h~'^ Mpc spline- 
smoothing kernel centered on the particle. We implement 
this with the SMOOTH code'^. From the distribution of 
the weighted particles, we calculate the density in each 
mesh in real space and redshift space. The finger-of-God 
effect is preserved in redshift space in this scheme. 

We first discuss the biased power spectra at z = 1 and 
z — 0.3 because the two share common bias schemes. We 
then present the z = 3 results. 

4.3. Bias effects at z = 1 

Figure 3 shows power spectra of biased tracers divided 
by a zero-baryon power spectrum in real space and redshift 
space. The left and right panels show number-weighted 
cases (NUM) and mass- weighted cases (MASS), respec- 
tively, with the minimum group multiplicity to, of 4, 10, 
and 30. 

NUM cases with to — 4, 10, and 30 produce biased 
tracers with h ^ 1.7, 2, and 2.5 where h is the ratio of 

''http:/ /www-hpcc. astro, washington.edu/tools/ 
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panel is for the number- weighted cases (NUM), and the right panel is for the mass- weighted cases (MASS). Green dashed line: the input 
power spectrum, black: the nonlinear matter power spectrum at 2 = 1, red: biased with m = 4, blue: biased with m = 10 and violet: biased 
with m = 30. 



the biased power spectra to the matter power spectra on 
large scales. The mass-weighted tracers (MASS) in the 
right panel of Figure 3 show larger biases than number- 
weighted cases for the same m: b ^ 2.4, 2.7, and 3.1. This 
is because the mass-weighted cases give more weight to 
the high-mass halos with larger biases while the number- 
weighted cases are dominated by halos close to the mass 
threshold. 

In both cases, the relative anomalous power increases 
as the amplitude of bias increases. For example, the 
m = 4 cases show little anomalous power, and the m = 10 
and m = 30 cases show larger but still mild anoma- 
lous power for k < /cmax (0.19/i Mpc^^). When different 
bias models with the same multiplicity are compared, the 
MASS cases have slightly larger anomalous power up to 
k ~ 0.3h Mpc~^. When different bias models with a simi- 
lar value of bias are compared (e.g., m = 4 of MASS and 
m = 30 of NUM) , the MASS case exhibits smaller anoma- 
lous power up to A: 0.3/iMpc~^ than the NUM case. 
MASS cases with all galaxies placed at the center of a halo 
showed the same trend, which means that the trend is not 
due to the effect of the halo density-profile. Thus, the bias 
scheme for the NUM cases introduces larger anomalous 
power, be it from the shot noise or from the nonlinear bias 
effect. One interesting question will be whether the differ- 
ence in anomalous power directly relates to the erasure of 
baryonic features. 

The upper panels in Figure 4 show d In P/d In k of biased 
power in real space (solid lines) in comparison to the input 
power spectrum (dashed fines). In the figure, all biased 
power spectra in real space preserve oscillatory features 
at least up to fc ~ 0.2ft, Mpc~^ while the contrast appears 
decreasing due to the nonlinear bias effect as bias increases. 
As before (§ 3.3), adding a smooth component due to a 
shot noise or nonlinear bias to the power spectrum will 
decrease the contrast in baryonic features. Subtracting off 



this smooth component will help to recover some of the 
baryonic features. We will revisit this in § 4.4. 

In redshift space (dashed lines in Figure 3), the power 
spectra of the MASS cases clearly show the finger-of-God 
suppression as k increases. The finger-of-God suppression 
is mostly removed in the spectra of the NUM cases, as 
is to be expected from the methods of biasing, and they 
show a mild remnant suppression in power with respect to 
the Kaiser formula^ as k increases. The lower panels of 
Figure 4 show that baryonic oscillations are smeared more 
in redshift space than in real space not only in MASS cases 
but also in NUM cases despite the suppressed finger-of- 
God effect in the latter case. Again, until we remove the 
bias effect on the broadband shape, it is hard to determine 
the degree of erasure. 

4.4. Bias effects at z — 1 with the broadband shape 
restored 

We next subtract the anomalous power from the biased 
power spectrum to restore the broadband shape and elim- 
inate the superficial decrease in contrast of the baryonic 
features. The restoring process we adopt is intended to 
assess the optimal amount of information on baryonic fea- 
tures available from the biased tracers. In real galaxy red- 
shift surveys, this unbiasing would be determined simul- 
taneously in the parameter estimation, which is an addi- 
tional complication. 

We fit the biased power spectra to a multiple of the non- 
linear matter power spectrum b^Pmattcr at the given red- 
shift plus a polynomial function /nl — cq+ci k+C2 k^. We 
then subtract the smooth function /nl that represents the 
anomalous power from the biased power spectra when cal- 
culating din P/d In A:. As the anomalous power from bias 

^Kaiser formula in our context means (1 -|- 2f3{k)/3 + /3^(A;)/5) 
where /3(fe) = Q^^{z)/b{k), and b{k) is calculated from the ratio 
between the biased power spectra and the matter power spectra. 
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Fig. 5.- din P/dln k from Pbiased — /nl(c(), fe, fe^) at z = 1 in real space. The left panel is for the number- weighted cases (NUM), and the 
right panel is for the mass-weighted cases (MASS). The label m denotes the minimum group multiplicity for selecting halos. Gray dashed 
line: the input power spectrum, gray solid line; the nonlinear matter power spectrum at z = 1, black solid lines: the biased power spectra at 
z = I. The vertical dashed line denotes fcmax = 0.19/i Mpc~^. One sees that the biased power spectra follow the features in the underlying 
matter power spectrum fairly well on linear and quasilinear scales once the broadband shape is restored. 



is smooth, this wiU help recover baryonic features. The 
ranges of wavenumber in real or redshift space is chosen 
suitably at each redshift so as to span well beyond the lin- 
ear region to constrain /nl but short enough so as not to 
weight the fit too much towards large wavenumbers. Vari- 
ations in /nl due to choosing different ranges of wavenum- 
bers or different degrees of the function (up to 3''''-order) 
do not show a meaningful impact on dhiP/dhik on large 
scales. The resulting function /nl behaves as constant on 
large scales as required, although on small scales, /nl for 
the MASS cases show a slow roll-over due to the extended 
halo profiles compared to the NUM cases. 

Figure 5 shows the resulting derivatives at ^; = 1 af- 
ter anomalous power /nl is subtracted (black lines) in 
comparison to the nonlinear (solid gray lines) and input 
matter power spectra (dashed gray lines). For k < fcmax, 
the agreement between the biased power spectra and the 
nonlinear matter power spectrum is excellent regardless of 
the different biasing schemes. Beyond fcmax, we observe 
small variations depending on bias schemes. The varia- 
tions appear related to the amount of anomalous power 
as this contributes additional power above the underlying 
baryonic features, increasing the statistical noise. Within 
the same bias models, higher mass thresholds and hence 
larger biases yield noisier derivatives, as would be expected 
from the increase in anomalous power. If we compare two 
bias models with similar bias values, m = 4 of MASS and 
TO = 30 of NUM, then the baryonic features in to = 30 of 
NUM appear noisier as it has a larger anomalous power. 
However, when those with the same group multiplicity are 
compared between different bias models, the baryonic fea- 
tures in the MASS cases look no worse than those in the 
NUM cases even with the larger bias and the slightly larger 
anomalous power for k < 0.3/i Mpc~^. 

Despite slight variations depending on bias schemes, the 
baryonic features in general have not been obviously dam- 



aged by the biasing process; the biased power spectrum 
closely follows the details in the underlying nonlinear mat- 
ter power spectrum over a broad range beyond /Cmax- This 
is different from the effect of nonlinear gravity shown in 
§ 3.3 where we saw the effect of mode coupling. This sup- 
ports the hypothesis that the anomalous power, whether 
due to shot noise or nonlinear bias, is a smooth function 
of wavenumber. The contrasts of the baryonic features 
appear slightly larger in some of the biased power spectra 
relative to the matter power spectrum, but this is likely 
due simply to the increased noise. 

For the redshift-space power spectrum of the MASS 
cases, we apply a similar fitting process to the one in § 3.3 
to restore the broadband shape both from the nonlinear 
bias and nonlinear redshift distortions. First we correct for 
the finger-of-God effect with a multiplicative function Ffog 
to match the biased power in redshift space to the biased 
power in real space. We then calculate and subtract /nl 
to remove the anomalous power. For the redshift-space 
power spectrum of the NUM cases, we do not correct for 
the finger-of-God suppression despite the slight deviation 
from the Kaiser formula. We calculate and subtract an 
additive /nl from the redshift-space power in this case. 

Figure 6 shows the resulting derivatives at z = 1 in 
redshift space (black lines) in comparison to the non- 
linear matter power spectrum in real space (gray lines). 
In redshift space, the contrast of the last feature before 
k = 0.2h Mpc~^ is smaller than in real space but still in 
good agreement. Beyond fcmax, we see the traces of bary- 
onic features although they look noisier than in real space. 
Again, the nonlinear effect of redshift distortions on bary- 
onic features in the NUM cases are no better than that in 
the MASS cases even though the virialized motions within 
the halos are suppressed in NUM cases, and this probably 
is related to the nonlinear effect on the velocity fields on 
large scales (Scoccimarro 2004). 
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Fig. 6. — d\nP/d\nk from Pbiascd — ,fNL(c(), k, k'^) at z = 1 in rcdshift space. Note that the MASS cases are corrected for the fingcr-of-God 
suppression beforehand. The left panel is for the number- weighted cases (NUM), and the right panel is for the mass- weighted cases (MASS). 
Gray; the nonlinear matter power spectrum at z = 1 in real space, black: the nonlinear matter power spectrum in redshift space (NO BIAS) 
and biased power spectra in redshift space. Note that the matter power spectrum in redshift space (NO BIAS) is also corrected with Ff^g to 
fit to the matter power in real space. The fitting range (fc < fef^t ) is indicated by the extent of the black lines. The rmmber-weighted cases are 
fitted to fcfit = 0.5/i Mpc"'^ but fcfit = 0.35/i Mpc~^ elsewhere. The vertical dashed line denotes the value of fcmax- The redshift-space biased 
power spectrum nearly reproduces the bajryonic features of the real-space biased power spectrum with mild degradation. 



4.5. Bias effects at z = 0.3 

Wc next investigate the effects of bias at lower redshift. 
Figure 7 shows power spectra of biased tracers divided 
by a zcro-baryon power spectrum at z = 0.3. The NUM 
cases (left panel) generate tracers with b ~ 1.2, 1.3, and 
1.6, and the MASS cases (right panel) generate 6 ~ 1.8, 2, 
and 2.3. At this redshift, we are particularly interested 
in tracers with 5^2, which corresponds to the lumi- 
nous red galaxy sample (LRG) of Sloan Digital Sky Survey 
(SDSS). This corresponds to m = 10 in the MASS cases 
although anomalous power in the m = 10 case is half the 
inverse of the number density of galaxies in the LRG sam- 
ple, suggesting that the m = 10 bias model is not exactly 
right. The bias and relative anomalous power is small up 
to A: ^ 0.15ft, Mpc~^ in all cases relative to ^ = 1. Recov- 
ering baryonic features beyond kjnax (O.ll/i Mpc~^) will 
be possible both in real space and in redshift space from 
Figure 8. 

Figure 9 shows the derivatives in real space after a 
corresponding smooth function /nl is subtracted. From 
the figure, the baryonic features in biased power spec- 
tra trace those in the matter power spectrum well up to 
k ~ 0.2/iMpc~^. Again, despite the larger biases of the 
MASS cases, they preserve baryonic features no worse than 
the NUM cases. Figure 10 shows the derivatives in red- 
shift space compared to the matter power spectrum in real 
space. All biased redshift-space power spectra trace the 
features in the real space matter power spectrum fairly 
well up to fc ~ 0.2h Mpc~^ but with more degradation in 
the NUM cases. 

To summarize, subtracting the smooth anomalous power 
helps to recover the contrast of the baryonic features at 
z = 0.3 as well. With small biases (6 < 2), the recov- 



ered contrast is comparable to the contrast in underlying 
matter power spectra even in quasilinear scales meaning 
that the nonlinear scales deduced in § 3.3 is valid despite 
biasing. Also the recovered contrasts do not seem very sen- 
sitive to moderate variations of biases, which is consistent 
with the results at z = 1. 

4.6. Bias effects at z = 3 

We next show the results at our highest redshift bin, 
z = 3. The number of simulation boxes used for this red- 
shift is 30, which is smaller than the other redshift bins. 
We generated three biased tracers: m = 4 for MASS, 
m = 4 for NUM, and Rho^ (Figure 11). The former two 
cases generate power spectra with b ~ 5.5 and 4.9, which 
are too high for Lyman break galaxies (Steidel et al. 1996). 
Correspondingly, the number density of these halos is very 
small, leading to significant noise in the power spectra. 
The Rho^ model, on the other hand, generates a bias of 
b ~ 2.5, similar to that of Lyman break galaxies (LBG). 
The anomalous power of Rho^ is 70% of the shot noise ef- 
fect from the number density of 10~^h^ Mpc~^ that we as- 
sumed for the sample in SE03 although the relative effect is 
nearly equivalent. The relative anomalous power amplifies 
the power by a factor of two at fc ~ fcmax(= 0.53/i Mpc~^) 
for the Rho^ case. 

While the matter power spectrum preserves baryonic 
peaks up to fc ^ QAh Mpc~^ (Figure 1), Figure 12 shows 
that even the Rho^ case cannot probe baryonic features 
beyond fc ^ 0.35/i Mpc~^ either in real space (left panel) 
or redshift space (right panel) . The MASS and NUM cases 
trace the acoustic oscillations only up to fc ^ 0.25h Mpc~^. 
That is, the tracers with very large bias do not mimic the 
underlying matter power spectra very well, unlike tracers 
with moderate biases at lower redshift bins. This is likely 
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Fig. 7. — Biased power spectra at z = 0.3 divided by a zero-baryon power spectrum. The left panel is for the number-weighted cases 
(NUM), and the right panel is for the mass- weighted cases (MASS). The label m denotes the minimum group multiplicity. For the NUM 
cases, b ~ 1.2, 1.3, and 1.6, and for the MASS cases, b ~ 1.8, 2, and 2.3 as m increases. The solid lines are for the real-space power, and the 
dashed lines are for the redshift-spaco power. 
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Fig. 8. — dlnP/dlnk from the biased power spectra at z = 0.3 in real space (upper panels) and redshift space (lower panels). The left panel 
is for the number- weighted cases (NUM), and the right panel is for the mass- weighted cases (MASS). Green dashed line: the input power 
spectrum, black: the nonlinear matter power spectrum at z = 0.3, red: biased with m = 4, blue: biased with m = 10 and violet: biased with 
m = 30. 
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Fig. 9. — dlnP/dlnA; from 

-^biased ^ /nl(coi ^i ^*^) at 2 — 0.3 ill real space. The left panel is for the number-weighted eases (NUM), and 
the right panel is for the mass- weighted cases (MASS). Gray dashed line: the input power spectrum, gray solid line: the nonlinear matter 
power spectrum at 2 = 0.3, black solid lines: the biased power spectra. Fitting to fegt = O.ShMpc"^. The vertical dashed line denotes 

fcmax = O.ll/l Mpc-^ 
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Fig. 10. — dlnP/dlnfc from i-^iased ~ fNh{co,k,k'^) at z = 0.3 in redshift space. Note that the MASS cases are corrected for the finger- 
of-God suppression beforehand. The left panel is for the number- weighted cases (NUM), and the right panel is for the mass- weighted cases 
(MASS) . Gray: the nonlinear matter power spectrum in real space, black: the nonlinear matter power spectrum in redshift space (NO BIAS) 
and biased power spectra in redshift space. The fitting range {k < fegt) is indicated by the extent of the black lines. For the number-weighted 
cases (NUM), fegt = 0.5h Mpc~^ but kfn = 0.35h Mpc~^ elsewhere. The vertical dashed line denotes the value of femax- 



12 



70 
40 



30 
10 



_l 1 1 1 1 1 1 1 1 1 


1 1 1 1 1 

MASS m=4 


1 1 1 1 i_ 




-0.5 










-1 


-1 1 1 1 1 1 1 1 1 1 


1 1 1 1 1 1 1 


1 1 1 1 1- 




-1.5 




NUM m = 4 




Mink 


-2 


1 1 1 1 1 1 1 1 1 1 


1 1 1 1 1 1 1 1 1 1 1 


1 1 1 1 1 


c 










■3 


-0.5 




Rho2 






-1 


M J^^l 1 1 1 1 1 1 1 


1 1 1 1 1 1 1 1 1 1 1 


1 1 1 1 1 




-1.5 


h No bic 


IS 






-2 


., , , 1 , ,V-n~rT 






0.1 0.2 


0.3 0.4 


0.5 







k (h Mpc" 



.Jl 1 1 1 1 1 1 

:\ Real 
! /A 

: ./ 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 L 

space ; 


= , , , , 1 , , , , V , , ,', 1 , , , , 1 , , , , 1 , ,- 



0.1 0.2 0.3 0.4 
k (h Mpc-i) 



0.5 



Fig. 11. — Biased power spectra at 2 = 3. Left: biased power divided by a zero-baryon power spectrum. Solid lines are for the real space 
clustering, and dashed lines with the same color are for the corresponding redshift-space clustering. Right: din P/dlnk from the biased power 
spectra. Green dashed line: the input power spectrum, black: the matter power spectrum at 2 = 3, red: Rho^ , blue: NUM with m = 4 and 
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because of the statistical noise from a small number of 
high mass halos, but we cannot exclude the possibility of 
an emergence of mode coupling as the nonlinear bias effect 
becomes very large. 

We estimate that galaxies with 5 ~ 3 at this redshift 
will recover baryonic features up to wavenumber of about 
0.3/i Mpc~^. The details of the result may vary depending 
on the biasing schemes and the shot noise. 

4.7. Summary of the effect of bias 

To summarize the general effect of bias on baryonic fea- 
tures, subtracting a smooth function to match the slope 
of the matter power spectrum largely restores the bary- 
onic features in the underlying matter power spectra if the 
amplitude of bias is moderate. This implies that a mod- 
erate anomalous power of bias whether from shot noise 
or a nonlinear bias effect does not erase the initial fea- 
tures. In other words, it is reasonable to assume that the 
anomalous power is smooth in wavenumber and does not 
generate features that mimic baryonic oscillations. 

We find that the detailed effects of bias not only scale 
with the amplitude of bias and the anomalous power but 
also depend on the biasing schemes used. For very large 
amplitudes of bias, we clearly lose information whether it 
is due to increased nonlinear bias effects or merely shot 
noise. 

Based on the baryonic features preserved in the matter 
power spectra (§ 3) and the effect of bias, we summa- 
rize the nonlinear scales in the biased power spectra. At 
z = 0.3, the results are encouraging in that moderately 
biased power spectra (6 < 2) preserve baryonic features at 
k < fcmax (O.ll/i Mpc^^) and record a fair amount of lin- 
ear information even beyond fcmax in both real and redshift 
space. The biased power spectra at z = 1 (6 < 2) show de- 
creased contrast for k < /cmax (0.19/i Mpc~^) and contain 
attenuated traces of baryonic features beyond /cmax- At 
z = 3, the biased power spectrum with 6 ~ 3 will preserve 



the features up to fc ~ 0.3h Mpc""'^, but no further. 

The biased power in redshift space traces the real-space 
biased power reasonably well with partial degradations in 
the baryonic features depending on biasing schemes. Sup- 
pressing the finger-of-God effect does not help to preserve 
real-space features any better, implying that the motions 
between halos do not strictly respect linear theory. 

To this point, we have shown the successful restoration 
of baryonic features in linear and quasilinear scales from 
various biasing schemes based on halo-mass thresholds. 
Given that all mass thresholds preserve the features, it 
is unlikely that more complex descriptions of galaxy pop- 
ulations, e.g., superpositions of our biasing schemes, would 
change the results. Similarly, a more stochastic local bias 
would not likely remove the oscillations, although stochas- 
tic models can increase anomalous power (Dekel & Lahav 
1999; Scherrer & Weinberg 1998) thereby reducing con- 
trast and increasing noise. 

5. IMPACTS ON COSMOLOGICAL DISTANCE ESTIMATION 
: ANALYSIS OF THE N-BODY DATA 

We next consider the impact of nonlinearity and bias on 
the statistical constraints on cosmological distances from 
the baryonic features. We are interested in using the power 
spectrum measurement to constrain a distance scale, and 
therefore we want to estimate how well we can constrain 
dilations in wavenumbers. We define this dilation parame- 
ter as a(= krcf / ktruc) ■ The error on the dilation parameter 
a represents the errors on the angular diameter distance 
Da{z) and Hubble parameter H{z). We perform a anal- 
ysis to fit the spherically averaged power spectra in real 
space from N-body simulations, Pobs, to a linear combi- 
nation of the input linear power spectrum, Puncar, and an 
additional polynomial function (eq. [1]). The fit parame- 
ters are a, a multiplicative bias bo, a scale-dependent bias 
bi, and additive terms for nonlinear growth or an addi- 
tional constant (oq, oi and 02). The mean value of a is 
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Fig. 12. d In Pj d In k from Pbiaaed — /nl(co, fc, fc'^) at z = 3 in real space and rcdshift space. The Rho^ case and the MASS case in redshift 
space are corrected for the finger-of-God suppression beforehand. Left: the real-space power spectra. Gray dashed line: the input power 
spectrum in real space, gray solid line: the matter power spectrum in real space, black solid lines: the biased power spectra in real space. 
Right: the redshift-spacc power spectra. Gray solid line: the matter power spectrum in real space, black solid lines : the biased power spectra 
in redshift space. Note that the y-axes in the left and the right panels are not scaled the same unlike the previous figures. The vertical line 
is at femax of 0.53/1 Mpc~^. From the figure, the tracers with a very large bias do not mimic the underlying matter power spectra very well. 
One sees that the mildly biased Rho^ case (6 ~ 2.5) probes baryonlc features up to fc ~ 0.35fe Mpc~^. 



expected to be unity since we set the reference cosmology 
to be equal to the true cosmology for simplicity. 

-fobs(^ref) — (^0"l"^l^ref ) -^-Plinear 

(1) 

We try two cases, with and without including bi. In 
both cases, we fit power at wavenumbers beyond fcmax to 
set the nonlinear trends. We note that even though the 
model spectrum contains higher harmonics of the baryon 
oscillations that may be absent from the data, this need 
not bias the distance measurement. The smooth portion 
of the nonlinear power spectrum that replaces higher har- 
monics does not have narrowband features to match on the 
linear model and hence impose a preferred physical scale. 
We will return to this point later in this section. 

The mean value and error of a are computed using 
jack-knife subsampling of the simulations (Lupton 1993). 
Since we do not know the true covariancc matrix, we as- 
sume a Gaussian error in each k-h'm constructed from the 
power spectrum averaged over all sets. Although assum- 
ing a Gaussian error implies that we underestimate the 
statistical noise relative to the true non-Gaussian error, 
the variations among jack-knife estimates of a in the 51 
subsamples (30 a.t z — 3) should reflect the non-Gaussian, 
mode-coupled error as these subsamples are drawn from 
actual nonlinear N-body data. In other words, our fit- 
ting slightly miswcights the data relative to optimum but 
should not produce overly optimistic aa compared to the 
true error. 

The behavior of resulting errors are consistent with the 
effect of nonlinearity and bias that we studied in the pre- 
vious sections. For the underlying matter power spectra 
at 2; = 3, we derive ac, ^ 0.35(0.36)% with wavenumbers 
less than fcfit = 0.3(0.5)/i Mpc~^ For the Rho^ biased 



power spectrum, ctq ^ 0.35% when k^t ^ 0.3/iMpc~^. 
Increasing the range of k beyond this value increases the 
error, which is consistent with the noisy feature in the left 
panel of Figure 12. If we scale the error to the survey 
volume assumed in SE03, this error value corresponds to 
1%. The anomalous power in the Rho^ case is close to the 
shot noise we assumed there. The analytic results in SE03 
implies aa ~ 0.93% (from 1/(7^ = l/o-|,^(^) + ^/(^%(^)), in 
good agreement. 

At z = 1, we calculate ctq ^ 0.4% for the underlying 
nonlinear matter power spectrum if we include a region 
up to fcfit = 0.3h Mpc~^. Among the biased power spec- 
trum, the NUM case with m = 4 has a bias value close to 
that of the redshift bin at z = 1 in SE03. For this case, 
(Ta ^ 0.4 — 0.5%, which corresponds to ctq ~ 1-6 — 1.8% 
when scaled to the baseline survey volume in SE03. This 
is to be compared with cr^ of 1.4% for the corresponding 
number density in SE03. Thus the simulation indicates 
a slightly worse precision relative to our previous predic- 
tions in SE03. The equivalent fcmax for a linear approxi- 
mation will be 0.17 — O.I8/1 Mpc~^ at 2; = 1, rather than 
0.19/?,Mpc-i. 

At z = 0.3, we find more optimistic results, as we would 
expect from Figure 9. For the underlying nonlinear matter 
power spectrum, we get CTq ~ 0.6% for kat = 0.3h Mpc~^. 
For the MASS cases with m = 10 and m = 30, which 
are similar to LRG samples, we find (Ta ^ 0.8 — 0.9%. 
When scaled for a survey volume of Ihr'^ Gpc'', CTq, is 2.1 — 
2.3% (the equivalent fcmax ^ 0.15 — 0.155/?, Mpc"-^), which 
is better than the 3.9% calculated from values in SE03. 
This is to be compared with the current observations: the 
4% measurements from Eiscnstcin et al. (2005) would give 
cTa ~ 3% when scaled to Ih"^ Gpc'^. The cause of the 
difference between 2.1 — 2.3% and 3% is due to the neglect 
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of redshift distortions in this modeling. 

In general, the errors calculated with and without bi are 
consistent. Without bi, the mean values of a are close to 
unity, indicating negligible bias. But a is slightly biased 
above 1 in cases with 61, particularly at lower redshift, 
albeit by < 1% in most of cases. Since we use the lin- 
ear power spectrum to match nonlinear power, the fitting 
process favors a negative 61 to match the erased baryonic 
features and, in order to compensate the resulting phase 
shifts of the oscillations, biases a preferentially above 1. 
On the other hand, without bi, the fitting process does 
not have means to use bias on a to account for the erased 
features, giving little bias on a. Including an appropriate 
recipe to account for the erasure of the baryonic features 
should remove this bias. That is, this bias would likely be 
easy to calibrate and remove using N-body simulations of 
a reasonable cosmological model. 

The results of analysis can be translated to the survey 
volumes required to achieve a a ~ !%• Table 1 presents 
estimates of the required survey volumes, assuming biased 
power spectra in real space. 

We do not extend the analysis to the redshift-space 
power spectra because of our lack of a reliable model of 
the redshift distortions to fit. There are deviations from 
the Kaiser formula on large scales as well as the finger- 
of-God effect on intermediate and small scales, both of 
which could involve an arbitrary angular dependence in 
two dimensions (reduced from three dimensions by the 
azimuthal symmetry). Nevertheless, the comparisons be- 
tween the real-space and redshift-space power spectra in 
our study lead to qualitative estimations of the effects of 
redshift distortions. Due to the decreased contrast of the 
baryonic features we observed in the spherically averaged 
power spectra in redshift space, we expect that the errors 
on H{z) will be degraded relative to the analytic predic- 
tion. At z = 3, we expect the effect will be insignificant. 
At z = 1, the degradation will produce a larger error on 
H{z) than in SE03. At z = 0.3, a degradation due to the 
redshift distortions will increase aH but likely no worse 
than the estimates in SE03. It is important to note that 
the information in the spherically averaged redshift-space 
power spectra will not be quantitatively equivalent to that 
in the two-dimensional redshift-space power spectra. Here 
we are averaging out the nonlinear redshift distortions, 
which are actually angle-dependent. For example, a sim- 
ple exercise of analysis using only the modes nearly 
along the line of sight suggests that the degradation of 
H{z) at z = 1 could be as large as a factor of two between 
real space and redshift space. We plan to investigate this 
in a future study. 

6. CONCLUSION 

We have used a large set of N-body simulations to show 
that the baryonic oscillations from the large galaxy red- 
shift surveys survive on large scales well despite the mild 
nonlinearity of gravity, redshift distortions, and bias. We 
compared the nonlinear effect on the baryonic features ob- 
served in the N-body results with the choices of nonlinear 
scale kynax in SE03. 

As expected, the nonlinear gravitational evolution erased 
the baryonic features progressively from smaller scales to 
larger scales as the redshift decreased. In real space, the 
nonlinear scales we have assumed in SE03 seem fairly con- 



servative at z = 0.3 and z = 1. We find that we need a 
slightly smaller fcmax for z = 3, but this modification has a 
minor effect on standard ruler test due to the small nonlin- 
ear scale. The redshift distortions imposed an additional 
obscuration for k < /cmax- Nevertheless, the redshift-space 
power spectra reasonably traced the features in the real- 
space power. 

We have shown that moderate nonlinear bias (6 < 3) 
does not erase the initial features. The biased power spec- 
tra follow the features in the underlying matter power 
spectrum fairly well once the broadband shape is restored. 
The effect of bias is not only proportional to the ampli- 
tude of the bias and anomalous power but also depends on 
the biasing models as well. However, these dependences 
seem to have minor effect on the underlying baryonic fea- 
tures. In redshift space, suppressing the finger-of-God ef- 
fect does not improve recovery of the features, and this 
indicates that nonlinear effects in the velocity fields on 
large scales obscure baryonic features as well. Neverthe- 
less, the redshift-space biased power spectrum reproduces 
the baryonic features of the real-space counterpart reason- 
ably well despite the degradation mainly due to the modes 
along the line of sight. 

From analysis of N-body results in real space, we 
predict errors on cosmological distances similar to those 
in SE03. Thus the nonlinear scale fcmax we have adopted 
in SE03, with minor modification detailed below, ade- 
quately describes the effect of nonlinearity on the standard 
ruler test. Furthermore, this implies that the cosmological 
distortions will be indeed distinguishable from nonlinear 
growth and scale-dependent bias, and so the derived un- 
certainty on cosmological distances depends on the degree 
of erasure of the baryonic features. Considering nonlinear 
gravity, redshift distortions and the clustering bias effect 
from N-body results all together, we estimate nonlinear 
scales appropriate for calculation of the information in the 
baryonic features. We consider both the loss of informa- 
tion for k < fcmax and the additional linear information 
on smaller scales which compensates the loss. Referring to 
the results from the analysis, the biased power spectra 
at z = 3 with b ^ 3 traces baryonic features well up to 
k ~ 0.3/1 Mpc~^ both in real and redshift space. At z = 1 
with b<2,k max ~ 0.17 — O.lSh Mpc ^ are appropriate in 
real space, but we need a slightly smaller /cmax for redshift 
space. At z = 0.3 with b < 2, k^^^^ ^ 0.15ft- Mpc"-^ for 
real space but smaller fcmax(> 0.11/;, Mpc^^) if we consider 
redshift distortions. We translate these results to the dis- 
tance measurements: at z = 3 and z = 0.3, we expect that 
Da{z) will be constrained as well as the estimates in SE03 
while H{z) will be slightly less well constrained because 
of the nonlinear redshift distortion effect on the baryonic 
features. At z = 1, the deviation will be the largest. While 
Da{z) will be near the estimates in SE03, H{z) can be as 
large as twice of Da{z). 

To summarize, the standard ruler test using baryonic 
features are robust against nonlinear effects in the lin- 
ear and quasilinear regime. Therefore, using the standard 
ruler, the on-going and future large galaxy redshift sur- 
veys will measure the dilations in observed scales due to 
cosmological distortions at various redshifts to excellent 
accuracy, providing a superb probe of the acceleration his- 
tory of the universe. 
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Table 1 

Required Survey Volume to Achieve aa = 1% 



Rcdshift 


bias 


V (h-^ 

•'survey V'*' 


Gpc^) 


Uesih^ Mpc-") 


0.3 


2.0 


4.4 




2.1 X lO""* 




2.3 


5.5 




1.2 X 10"^ 


1 


1.7 


1.4 




2.1 X 10-3 




2.4 


1.7 




3.4 X 10"^ 


3 


2.5 


0.5 




1.4 X 10-3 



Note. - The approximate survey volumes required to achieve <7a = 1% for the biased power spectra. Note that the values are based on our ; 
analysis of the real space power spectra, and this means that redshift distortions are not accounted for. The effective number density rieff 
the inverse of the anomalous power at fc ~ 0. 



